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Comment on: "Pipe Network Model for Scaling 
of Dynamic Interfaces in Porous Media" 

In a recent Letter Q , Lam and Horvath (LH) present 
a pipe network model of imbibition. The network was 
made of small pipes uniformly distributed in sizes and 
large pipes of a single radius. By tuning the size and 
relative number of large to small pipes LH were able to 
find power law exponents close to those measured in the 
experiment of one of the authors (2| . For example model 
parameters could be chosen such that the mean velocity 
of the front (v) is related to the mean height of the front 
(H) as v ~ H~ n with Q = 1.6. A new exponent identity 
is proposed (i.e., Eq. (7) of the Letter) relating tempo- 
ral and spatial correlation functions. In this comment 
we point out that LH neglect an important length scale 
and that their exponent identity, whose underlying as- 
sumption holds trivially for standard kinetic roughening, 
should not apply to imbibition. 

To explain the main flaw in LH's arguments we recall 
the picture of spontaneous imbibition developed in Refs. 
HQ . Competition between the effective line tension and 
liquid transport imposes to the front a lateral length scale 
£ x regardless of the value of f2 (indeed, in experiments 
using paper as the porous material a variety of fi's dif- 
ferent from unity and a variety of scaling behaviors have 
been reported H). £ x is important since for any fixed 
height the width of the interface is finite and does not 
obey a power law relationship with system size or time. 
The width only diverges if the average height of the in- 
terface diverges as is the case for freely rising fronts. For 
example in the model of Ref. j|] the width follows w ~ 
with a global roughness exponent x — 1-25 and the lateral 
correlation length scales as £ x ~ v~ x l 2 ~ i 1 / 4 . 

The existence of an intrinsic lateral correlation length 
becomes obvious when systems of different sizes are con- 
sidered. Although such an analysis was not presented by 
LH a lateral correlation length can be inferred from their 
description of the height difference correlation function 
given as; 


C(l) = v- K g(lv iei+K ^ a ), 


(1) 


where g(x) ~ x a for x -C 1, and constant for large x. The 
reported values of the exponents are k = 0.49, a = 0.61, 
and 9 1 = —0.25. Equation (|j) clearly defines a length 
scale ~ « -7 , where 7 = {61 + n)/a ~ 0.4, which then 
implies C{1) ~ Qg(r/£, v )- Thus, ~ v~ 0A is analogous 
to the £x ~ v^ 1 / 2 in Ref. S. Equation (El) also defines 
a global roughness exponent, \ — + = 1.25, 

which implies a "superrough" interface with anomalous 
scaling 0] , as could be verified from the structure factor 
S(k,t)=_(\h k (t)\ 2 ). 

To derive their exponent identity, LH assume that the 
dynamics are controlled by a single time scale which can 
be obtained by two means. For a moving interface, a 


width w is reached after a time t\ ~ iy( sl + 1 )/( KO ). On the 
other hand, another time scale can be obtained from the 
interfacial fluctuations when the average interface height 
above a reservoir is kept fixed. In the steady state, LH 
write the time correlation function in the form 


C(t) 


7(*« ( 


(2) 


with the scaling function fix) ~ x^ for x <C 1 and 

constant at large x. This suggests a time scale ti ~ 
v -(8 t+K )/l3 _ W (e t+K )/ K p^ where Q t = 38 and p = 63 

are found. Assuming that t\ cx t%, LH obtain — 
tt(O t + k)/(Q + 1) » 0.54 which differs from the direct 
numerical fit quoted above by ~ 14%. 

Our main objection to the new exponent identity is 
that it lacks justification since t\ and ti describe two 
distinct physical time scales. It takes a time t\ for the 
correlation length to have value £, v {ti)- This in turns con- 
trols the width as w ~ £*. On the other hand, t% is the 
relaxation time of the fluctuations when the interface is 
kept at a fixed height. In this case, the correlation length 
£„ is a predetermined constant and the saturation within 
this zone is obtained when the spatial extent of the cor- 
relations equals £„. The scenario proposed by LH cannot 
occur if ti <C t\ and indeed their data indicate ~ w 2 8 
and t\ ~ ui 3 ' 25 which shows their assumption to be false 
asymptotically. Notice that the exponents are close, so a 
reasonable range of scaling is needed to see the difference 
clearly. 
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